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ABSTRACT 


bPe-opeCUEdie scl ellite thresholds Eor precipitation spec- 
ification are explored with visual and infrared satellite 
data collocated with Service-A hourly observations for 137 
surface stations in the southeastern United States. The 
data span the month of August 1979 and total 70,623 observa- 
tions, including 538 daylight precipitation observations. 

The distributional and statistical differences of four 
satellite resolution sizes ranging from 484 to 2025 nmi2 are 
explored and determined to be significant in the representa- 
‘Jon of weather conditions. Precipitation and no-precipita- 
“ion data can be statistically differentiated with the 
visual and infrared mean and standard daviation values. 

For overcast ceiling reports, a Simpl linear bi-spec- 
tral threshold based on a 50% probability of precipitation 
is defined as extending from albedo 1.00 to 0.60 with asso- 
Clated cloud tcp temperatures 290K and 210K, respectively. 
For overcast and broken ceiling reports, an albedo greater 


than 0.80 specifies a 50% probability of precipitation. 
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I. INTRODUCTION 


The determination On pre Cipiltation occurrence and 
amounts is an important factor in scientific, commercial, 
and operational endeavors. Scientific uses are concentrated 
in the fields cf meteorolcgy, hydrology, and oceanography, 
where precipitation is essential in analysis, diagnosis, 
prediction, and verification. Within meteorology, precipi- 
tation serves as both a forcing and response ¢lement in the 
study of daily weather and climatology. Indeed, precipita- 
tion is a critical input for climate reasearch and into gen- 
eral circulation nodeéels which premise to extend the time 
frame of skiliful weather forecasts. Commercial uses encon- 
Peres agricuiture, forestry, transportation, communicaticns, 
Water resource management, and many others. 

Despite the importance of precipitation data to a vari- 
ety of fields, there are serious shortcomings in current 
precipitation determination. These shortcomings are due to 
areal and e¢conomic limitations imposed upon <+he land-based 
rainfall monitoring systems. A possible solution is embod- 
ied in precipitation information extracted from satesllicxe 
Gidtea . With the advent cf high resolution, multi-spectral 


Channel satellites inthe late 1970's, satellite derived 
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precipitation data are being studied as a viable method <o 
complement and supplement conventional rainfall data. 

The satellite image interpreter does 2 subjective analy- 
sis based on the gray shade variations, répresenting @ range 
of digital counts, that appear in the satellite image. How- 
ever, satellite data contain more information within the 
digital values than can be resolved by the human eye in pho- 
tographic images. The satellite digital counts input into a 
computer allow use of the full range of the digital values. 

Until recently, computer processing of satellite data 
has been confined largely to research uses. Acquisition, 
Storage, and precessing of the huge volumes of digital sat- 
ellite data could not be handled operationally in real time. 
However, with the recent advent of more capable mini-con- 
puter systems, such as the United States Navy's Satellite 
Data Processing And Display System (SPADS) developed by the 
Naval Envircnmental Prediction Research Fracility (NEPRF) 
Monterey, California, real time quantitative use of digital 
Satellite data has become @ reality. With the operational 
availability of such systems as the SPADS unit, there is a 
need for numerical schemes to aid in the objective 


Specification of current weather conditions. 
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This thesis concentrates on the specification of visual 
and infrared satellite data thresholds in det¢rmining pre- 
Smpptacion occurrence and qualitative precipitation rates in 
a mid-latitude coastal environment. Thea data set used con- 
sists of collocated Geostationary Operational Environmental 
Satellite--East (GOES-E) satellite data and hourly surface 
observations at East Coast and Gulf Coast United States sta- 
tions, south of 40°N, for the month of August 1979. 

The use cf satellite data fer precipitation sp¢ecifica- 
—pon 1S not new. There is the recognized limitation that 
infrared and visual satellite sensors ars measuring proper- 
t1les associated with eres etega particles and not precipi- 
meton sized particles. Nonetheless, Muench and Keegan 
me 9) specified quantitative precipitation rates, Liljas 
(198i1a, 1981b) specified qualitative precipitation rates, 
and Lovejoy and Austin (1579) delineated rain versus no-rain 
cases uSing visual and infrared satellite data. Del Beato 
(1981) used cloud top temperatures within a restricted cloud 
case classification to derive qualitative precipitation 
rates. 

The currently available precipitation study results are 


based on data sets with region, season, and size 
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Pema tacions. This reséarch effort will use data fren 
stations covering more than 420,000 sguarS nautical miles 
(nmi2) in the eastern and central United States with a total 
ef 70,623 observations (538 precipitation observations). In 
comparison, the relatively comprehensivs precipitation study 
of Muench and Keegan (1979) was based on 552 cases (300 
tainfall cases) from five stations in the northeastern 
United States for April through November 1977. The signifi- 
cantly larger size of the present sample will allow better 
statistical determination of apprepriate distributions and 
threshold value significance. 

The primary objective of this thesis 1s to investigate 
Beecitication of precipitation versus noO-precipitation from 
satellite visible and infrared digital counts. Addition- 
ally, in precipitation cases, the feasibility of qualitative 
specification cf light versus moderate/heavy precipitation 
and guantiative specification of convective versus continu- 
ouS precipitaticn are investigated. 

The thesis is organized into five chapters. Chapter II 
reviews the satellite data based pracipitation studies. 
Chapter III describes the data set, th= data processing and 


the testing progran. Chapter IV describes the results. 
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Chapter V states the conclusions and suggests further 


research. 
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A. INTRODUCTION 

The specification of precipitation and the estimation of 
rainfall rates using satellite imagery have been studied 
using a variety of methods over awide spectrum of time 
scales. This review will concéntrat? on thos¢ methods 
developed for synoptic scale and mesoscale analysis of pre- 
cipitation on a diurnal or shorter time scale. The methods 
reviewed include bi-spectral and infrarei threshold (Muench 
and Keegan, 1979; Liljas, 198ta, 1981tb:; Lovejoy and Aus- 
rn, 1979; Del Beato, 1981; Wylie, 1982) and life history 
Mecotield, 1981; GElttichect al, 19753); ScOme setae lai, 


fee os Wylie, 1979; Negri and Adler, 1981). 


B. BI-SPECTRAL AND INFRARED THRESHOLD 

The bi-spectral threshcld method, in which infrared and 
visual satellite data eee cea: involves mapping the extent 
and distribution of precipitation. Combining the visual and 
infrared data provides infcrmation on the cloud temperatures 


(infrared data) and on the cloud thickness (visual data). 


Thus, while use of the visual or infrared data alone may 
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have limitations in specifying precipitation, the combina- 
tion of both may succeed at specifying precipitation. The 
multi-spectral satellite channels introduced on satellites 
in the late 1970's yielded the possibility of bi-spectral 
thresholds. Threshold values and study condition parameters 


of selected bi-spectral studies are summarized in Table I. 


TABLE I 


Summary of Bi-spectral and Infrared [Threshold Values 


TIME OF THRESHOLD VALUE 
STUDY LOCATION CASES YEAR INFRARED VISUAL 
Muench and Northeastern 552 April- 23°C 0.60 
Keegan United States obser- November 
(1979) vations 1977 
Liljas Scandinavia = May 1979 =12 76 ito = 
(1981) August 15°C 
1979 
Lovejoy and Montreal 17 days June 1977 Slee. 80, 83" 
Austin =76-6; 
(1979) 41°C 


* Visual threshold based on normalized scale from 0 - 1 


Muench and Keegan (1979) studied precipitation specifi- 
Cation using GOES visual andinfrared satellite data and 
hourly rainfall climatological data for five stations in the 
northeastern United States for the period April through 


November 1977. Phe iF data set consisted of SoZ 
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observations, comprised of 300 rainfall observations and 252 
cases of either nonprecipitating cloudy or fair weather 
observations. The visual (1 km resolution) and infrared (7 
km resolution) satellite data were area averaged over 7 x 7 
square kilometers (km?) and 14 x 14 km2, respectively. A 65 
point visual data array (8 x 8 plus the center point) anda 
17 peint infrared data array (4 x 4 plus the center point) 
were centered cver each station. The GOES visual data were 
normalized using reflecticn values from Liou (1976) with the 
modification of lower absorption and higher transmission to 
compensate for Liouts treatment of the complete solar spec- 
trum. The anisotropic radiation of clouds was corrected 
With functions calculated by Muench and Ke2gan (1979) from 
ground-based radiometers and satellite neasurenencts. From 
these data, they determined probabilities for precipitation 
greater than .01 and .10 inches in one hour and the amount 
See precipitation for the hour following the satellite 
observation (see Figs. 1, 2, and 3). 

Muench and Keegan (1979) did not provide th2 standard 
deviations for the data in these figuras. However, they 
Stated there was "considerable uncertainty in the specifica- 


tion of rainfall amount." As an examples, they stated <hat 
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Figure 3. 


feeng Fig. 3 "for a one~hcur rainfall specification of 0.10, 
two-thirds of the values would fall between 0.25 and 0.04." 
Muench and Keegan stated that their figures emphasize the 
requirement fcr both visual and infrared data to specify 
precipitation amounts. 

Liljas (1981a, 1981b) developed a bi-spectral cloud 
Classification based on visual and infrared data from the 
polar orbiting TIROS--6 satellite (see Fig. 4 and Table II). 
The data set ccnsisted of a limited number of daily observa- 


tions, chosen for their synoptic characteristics, in May and 
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D ; HO 
VISUAL 
o Dimension Decision Space for Typing Clouds 
m Visual and Infrared Digital Counts (Table II 


a the symbols for the clouds) (from Liljas, 
a 


wm Ch 1 


Figure 4, 


aro 


August 1979 over a region encompassing Norway, Sweden, fFin- 
land, and the Baltic Sea with weather charts providing the 
feound truth. Based upon the precipitation threshold 


results of Muench and Keegan (1979), Liljas chose a cloud 
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TABLE II 
Gholawebacsweleatton £6 be used with Fig. 4 


(from Liljas, 1981a) 


Main and Cloud Types: 


1. 


2. 


3. 


4. 


5. 


6. 


7. 


8. 


Cumulonimbus 


Himbostratus 


Cirrostratus 


Cumulus congestus 


Stratocunulus 


Haze/Stratus 


Water 


a 
= 
_ 
b 
b 
l 
by 
¢ 


we 


fe 


d 
d 


~ 
Nw er 


@ mr 


NH 


£ 
£ 


4 


a, 


f; 


b 


Storm cloud with high top 
Squall cloud with scattered showers 


Large vertical thickness 
Rather low topside 


Dense cirrostratus 

Cirrus 

Thin cirrus over water 
Dense altostratus 

Thin altostratus over water 
Thin altostratus 


Dense altocumulus 
Large piled up cumulus 
Rather small and flat c:-mulus 


Dense stratocumulus 
Ordinary 

Slightly piled up cumulus 
with clear areas in between 
Very dense haze/stratus 
Dense haze/stratus 


Ordinary haze/scratus 
Cumulis humilis 


Haze over water 
Planting season spring or autumon 
Warm green season 


old 


= © 
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moo temperature threshold of -129C to -15°9C to classify 
cumulonintus and nimbostratus clouds. Sittacting., wlth this 
cloud classification and the assumption that the highest and 
densest clouds produce the maximum precipitation amount, 
Liljas suggested a qualitative precipitation intensity scale 
based on the sum of the visual and infrared satellite digi- 
*al counts (see Table III). Thes¢ sums represent the areas 
of the Liljas nimbostratus and cumuloninbus cloud types in 


his bi-spectral cloud classification (s3223 Fig. 5). 


TABLE Iifl 


Threshold Values Describing Precipitation Intensity Levels 


as Applied in Fig. 5 (from Lilj3is, 1981a) 


The Sum of Digital Levels 


Ch 1+ Ch 4 S93 Osc ligh = rain 
311-330 
331-350 
35 1-370 
37 1-390 
390 very strong rain 


Levejoy and Austin (1979) studied rain mapping of cloud 


areas based on GOES visual and infrared satellite data over 
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VISUAL 


Precipitation Intensity Classification from 

Visual and Infrared Digital Counts. The 
Precipitation Area is Represented with the Dark 
Diagonal Lines. (See Table III for the 
mathematical description of the intensity areas.) 
(from Liljas, 1981a 


Canada, and the tropical Atlantic (Global Atmos- 


pheric Research Program Atlantic Tropical Experiment, GATE, 


data) with radar data providing the ground truth. The Mont- 


real data set consisted of 17 observations over three days 
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during June 1977. Working with 4 x 4 km resolution sa<el- 
lite image, Lovejoy and Austin plotted two dimensional fre- 
quency grids for the radar-determined rain and no-rain 
points on a 25 x 25 array (see Figs. 6 and 7). The visual 
data were normalized by selecting the "brightest" and "din- 
mest" values in each image and linearly interpolating the 
radiances kéetween 0 and 1. 

Lovejoy and Austin (1979) state, with reference to the 
mma@lus Frain data distribution of Fig. 6 that, "The distri- 
bution was to a good approximation a two-dimensional Gaus- 
sian." They do not describe or provide the statistics to 
Support this assertion. The no-rain cumulus cases (Fig. 7) 
were described as a bimodal distribution with one peak near 
the low visual and low infrared values and the other peak 
near the rain peak but shifted slightly toward lower values. 
In most cases, the separation of the cumulus rain and no- 
meee caseS WaS Statistically significant with the probabil- 
4ty ranging from 10% to 50% that che rain and no-rain 
Samples came from the same population. 

The Lovejoy and Austin (1979) two dimensional frequency 
plots for non-cumulus Storms were limited to one case. The 


Significant differences between the cumulus and non-cumulus 
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data sets were that the non-cumulus no-rain plot lost its 
bimodal character, relative to the cumulus no-rain plot, and 
appeared as a broad two dimensional Gaussian distribution. 
The non-cumulus rain plot points fell within the no-rain 
distributicn, but were shifted slightly higher in the vis- 
ual. The separation of the non-cumulus rain and no-rain 
cases was not statistically significant, with greater than a 
50% probability of the rain and no-rain samples coming from 
the same population. 

Lovejoy and Austin (1$79) attempted to further classify 
the cumulus rain and no-rain cases into no-rain, light rain, 
and heavy rain. Rainfall rates greater than 2 mm-h-!, as 
determined by radar, were defined as heavy train. As 
ex pected, the méan of the heavy rain cases was shifted 
slightly towards higher visual and infrared values than the 
mean of the light rain cases. However, the shift was so 
small that there was at least an 80% probabiiity of the 
light rain and heavy rain cases coming from the same popula- 
eavON Lovejoy and Austin (1979) concluded that "little if 
any rainfall-rate information is contained in a single (vis- 


ual and infrared) satellite image." 
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Lovejoy and Austin (1979) tested a2 spectral threshold 
technique for rain area mapping. Each satellit2 image of 
400 x 400 km was divided into one hundred 40 x 40 km boxes. 
The 100 sukareas were each checked with radar to determine 
the total number sf rain areas. An equal total number of 
satellite subareas were classified as rain areas. The sat- 
ellite subareas with the highest visual and highest (cold) 
infrared values were classified as rain areas, until the 
total number cf satellite rain areas equaled the total nun- 
ber of radar determined rain areas. This spectral threshold 
technique was applied to three days accumulation of data and 
is shown in Tables IV and V. When compared with «he success 
of the two dimensional frequency plot method, the visible 
and infrared thresholds averaged 45% and 58% worse, respec- 
tively. The accuracy of the visual threshold is limited by 
the extent of low, thick clouds and the infrared threshold 
is limited by the extent of the cirrus clouds in the satel- 
lite image. Lovejoy and Austin (1979) concluded that "the 
errors involved in using a 'best threshold' are very large 
indeed." 

Del Beato (1981) studied correlations between cloud top 


temperatures (based on NOAA-5 satellite data) and rainfall 
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TABLE IV 


Statistical Comparison of Rain Area Mapping Techniques 


(R /R X 100 indicates "percentage of correct satellite rain") 
(from Lovejoy and Austin, 1979) 


= 











Opt. 2-D Boundary 1R Optimum Threshold Visible Optimum Threshold 
ee a Area Rain Total No, 

Day (Ra/R) x 100° IR(K) (a/R) x 100 (Scale: 0-1) (Raj R):* 100 Coverage (%) of Points, 
a i Sess i i SS gS i SS 
GATE 
242, 243, 246 65 < 232 $3 > 0.68 : os 15.8 47706 
247, 248, 258 
252, 261 
Monrreal 
1S2 536 <232 20 >0.88 i 9.7 . 40361 
190 56 < 247 88 >0.80 aa 24.0 33738 
153° 33 < 234 $2 >0.88 ay 15.9 22538 

TABLE V 


Statistical Comparison of the Accuracy of Rain Areas 


(from Levejcy and Austin, 1979) 


a aaa 





Number of Images Error 
Technique Region or Sequences Bias Factor cams 
2-D Pattern Montreal 17 1.13 1.26 0.22 
Matching 
2-D Pattern Montreal 3 1.08 1.19 0.18 
Matching 
Optimum IR Montreal : 3 1.38 1.74 0.71 
Threshold ; 
Optimum Visible Montreal 3 1.54 1.59 0.58 
Threshold 
2-D Pattern GATE 8 Zi 1.4] 0.25 
Matching 





totals for 30- and 60-min intervals over eastern Australia. 
The satellite data had a 60 km2@ maximum resolution at subsa- 
tellite pcint and cloud top temperatures were area averaged 
for a resolucicn of 200 kme. The 21 data sets were first 


Meassified according to synoptic situation in a rough 





ae empe c©O gGEOUp wthe data by cleud typ2, droplet spectra, 
and air mass trajectory. The initial results suggested that 
+he cloud top temperature determined an upper limit on rain- 
fall amount, with the maximum increasing as the cloud top 
temperature decreases. A linear correlation analysis to 
determine a quantitative relationship between rainfall 
amount and cloud top temperature gave indefinite results. 
Further study of surface and radissonde observations 
indicated that classification by proportion of cumuliform 
cloud reports to all cloud reports and subcloud layer humid- 
ity might be mcre appropriate (Del Beato, 1981). This clas- 
Bemereation resulted ina correlation cosfficient of 0.90, 
excluding cases with cumuliform portions less than 50% and 
dew-point depressions of greater than 69°C. Finally, a con- 
posite frequency distribution was calculated based on three 
cases, all scuthwesterly stream situations described as 


Meest-frontal cellular convection cases in cyclonically 


curved flow." The fitted equation was: 


f = 0.057 - 0.004CTT - 0.054R (1) 


where f is the rainfall frequency, R is the 30-min rain 


eecal (mm), and CTT is the cloud top temperature (°C). The 


32 





equation was fitted to 41 independent f values. This equa- 
tion is associated with a correlation coefficient of 0.79 at 
+he 99% confidence level. Equation (1) indmweates RO rain 
from clouds warmer than +139°C and a maximum 30-min rainfall 
meee.) mm for a cloud top temperature of -20°C. 

In summary, Del Beato (1981) found that cloud top temp- 
eratures and 30- and 60-min rainfall totals indicated sta- 
tistically significant relationships for cloud systems with 
mmrign PLOPCrtion of cumulus clouds and high subcloud humid- 
ity. Additionally, as cloud top temperatures decrease to at 
least -35°9C, rainfall totals increase. 

Wylie (1982) attempted to correlats rainfall occurrence 
With radiosonde soundings, hourly Service-A observations, 
od visual and infrared satellite data. His data sample was 
restricted to "large-scale cloud cover" areas with wide- 
spread precipitation (rain gauge reports varied less than 
20%) for the Great Plains States region for the period 27 
February 1981 through 4 January 1982. From thirteen parame- 
ters derived from the three data sources (s2e Table VI), the 
best linear regression equation for e2stimating rainfall 


Tates was: 


6-hour rain (in) = 1.0242 + 0.380Pw - 0.03040c 


- 0.0047Ct (2) 
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where Pw is the vertically integrated pracipitablé water 
vapor (in), Qc is the moisture convergence (g/kg/day), and 
ct is the cloud top temperature (Kelvins). Equation (2) has 
a linear correlation coefficient of 0.60. Linear regression 
equations were also determined for the three parameters 
alone and for a combinaticn of Pw and Qc to be used when not 
all three data types were available. The cloud temperature 


regression eguation was: 


6 hour rain (in) = 2.10 - 0.008Ct (3) 


The correlation coefficient was -0.35. Wylie (1982) stated 
that the synoptic scale data base méasureaments were best 
suited for estimating broad changes in r2infall rates asso- 
Ciated with changes in air masses and not suited for esti- 
Mating rainfall rates asscciated with small scale dynamic 


processes. 


Ce mers HISTORY 

The life history methcds are empirically derived precip- 
itation estimation schemes based upon two assumptions, 
mies t , that significant rainfall comes from convective 
clouds, and second, that convective clouds can be identified 
and measured in satellite images. These methods involve 


manual analyses of convective cloud aréas in a sequence of 
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TABLE VI 


jmomemetatich Coefficients fer Determination of Precipitation 


Based on the Three Data Types (from Wylie, 1982) 


CORR ELATION 

PARANETS & PRECIP. REPORT OF CASES 

Vertically integrated 0.48 196 5 8* 
precipitable water vapor 

Cloud top brightness -0.44 184 4 yx 
Cloud top height -0.40 190 36* 
Moisture convergence Uns 184 31% 
Cloud top temperature =O a) 199 27* 
Bubble model predicted cond. Kaz! 115 9 

500 mb vorticity advection =) 272.1 N73 8 * 
Parcel lifted index =ia20 200 8* 
700 mb temperature advection Or 20 173 7* 
Sfc temperature advection Oe 156 6 

850 mb temperature advection Oey 189 6 

Wind convergence (sfc) 0.09 167 1 

Vertical wind shear O03 156 0 


Mao~tagniticant ccrtrelation at the 99% Lavel. 


visual, infrared, or both visual and infrared satellite 


images. Threshcld values and study condition parameters 
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associated with published life history studies are summa- 


rized in Table VII. 


TABLE Vil 
Summary of Life History Threshold Values 


TIME OF THRESHOLD VALUE 
STUDY LOCATION CASES YEAR INFRARED VISUAL 
Griffith, Florida, 34 days summers =20<6 80 counts* 
et al., Venezuela, 1969-1976 
(1978) Honduras, 
and hurricanes 
impacting East 
Coast United 
States 
Stout, tropical North 57 obser- September 226°C 0.45 albedo 
et al., Atlantic vations 1974 (sun over- 
(1979) head) 
Wylie Montreal 6 days June 1977 =16-G = 
(1979) September 
1977 
Negri and Oklahoma, 1 day April 24, a7 = 
Adler Arkansas, (15 thunder- 1975 
(1981) Missouri storms) 


® ATS-3 satellite 


The ScofieldyOliver (Scofield, 1981) analysis follows a 
decision tree procedure to estimate half-hourly rainfall for 
deep convective systems within tropical air masses. Using 


enhanced infrared and high resolution visual satellits 
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images, the technique involves first identifying the active 
Senvective portion Gf the cloud, or cluster, from two con- 
secutive satellite images. Once the active portion is iden- 
ca tied, the half-hourly rainfall estimation is computed 
based on such factors as cloud top temperature, G2oud 
growth, and departure of precipitable water from a summer- 
time normal. 

The Griffith/Woodley (Griffith et al., 1978) technique 
is designed =o estimate rainfall in the tropics, over large 
space and time scales, using geosynchronous visual or infra- 
red satellite imagery. This time-dependent technique was 
empirically derived as a relationship between cloud area, 
echo area, and rain rate for two aréas in south Florida, 
With raingag¢-radar providing the ground truth, and was then 
Bested in other tropical areas. This scheme was subse- 
quently tested further in extratropical areas (Griffith et 
als» 1980), with modifications to th¢ rainfall amount 
predicted. 

The determination of a cloud area-rainfall relationship 
first required the specification of both a visual and an 
infrared threshold to define the cloud area. The visual 


brightness threshold, normalized for radiation geometry, was 


of 





pom colmes fOr =hnep third Application Technology Satellite 
(ATS-3) and the infrared threshold was 253K (-20°9C). The 
thresholds were based on a comparison of the clouds with a 
given maximum digital count and the radar echoes associated 
with these clcuds. 

The empirical cloud reav-rainfall relationship was 
derived as a two st¢p process. Etictye a LelacLonship 
between the cloud area and the radar echo area, normalized 
for the maximum area achieved by the cloud or cluster, was 
established fer the visible and infrared satellite data. 
Second, the relationship between the echo area and rain vol- 


ume was determined and was of the form: 


Rv = I Ae (5) 


where Rv is rain volume per hour (m3-h-!), I is rain in units 
of (m3-km~*-h-!), and Ae is the echo area (km2) defined by 
the 1 mm-h-! rain rate. Thus, given a cime sequence of con- 
vective clcuds (or cluster areas) measured from visible or 
infrared satellit? images, volumetric rain rate can be 
estimated. 

Stout et al. (1979) mcdified the Griffith/Woodley tech- 


nique (Griffith st al., 1978) to estimate volumetric rain 
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rate directly from a cumulonimbus cloud area and area change 


according to the equation: 
R = aA + a,dA/dt (6) 


where R is the volumetric rainfall of the cloud (m3-s7-!1), A 
is the cloud area (m2), dA/dt is the change of cloud area 
over time (m@-s~1!), and a, and a,are constants with dimen- 
sions m-s-! and m, respectively. The two constants were 
calculated by a least squares fit of cloud area-rain rate 
pairs based on visible and infrared geosynchronous satellite 
@eea and 5.3 cm ship radar rain data collected during GATE. 
The cloud area and its change are defined by the threshold 
value. The visible threshold for cloud arsa calculations 
Was 60 digital counts on the ATS-3 (corresponding to an 
albedo of 0.45 with the sun overhead), or equivalently 172 
digital counts onthe first Geosynchronous Meteorological 
cece llite (SMS 1). The infrared threshold was 160 digital 
mouncs (~269C). The standard error betwe2n the estimated 
rainfall and tke mean radar rainfall was 62% and 76% for the 
visual and infrared equations respectively. 

Wylie (1979) attempted to use the tropical convective 


@eeetall techniques of Griffith et al., (1978) and Stout et 
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me, (tou) for estimating precipitation ian Montreal, Can- 
ada. Using visual satellite data, corrected for the chang- 
ing sun angle (Mcsher, 1975), infrared satellite data, and 
10.0 cm radar meaured rainfall rates, Wylie studied six days 
mee Drecipitaticn, three days each in June and September 
wa77 . Wylie concluded that because of air mass differences 
between Montreal and the tropics, ChS Greaetichs and Stout 
estimation techniques did poorly in Montreal, Canada. The 
Singlemost important limitation with these two schemes was 
the difficulty of measuring cumulonimbus cloud area when the 
"anvils were often merged into large cloud masses and the 
extensive stratus cloud cover often obscured the pictures." 
Wylie also noted that the Griffith el. al. (1978) threshold 
of -26°C had to ke changed to -16°C for the summertime Mont- 
real, Canada, area. With the warmer cloud top temperatures 
the cloud areas were a larger, more appropriate size for 
ma acking. 

Wylie (1979) then attempted to combine sounding data 
input inte a one-dimensional model (Simpson and Wiggert, 
1969) and satellite cloud cover measurements to estimate 
rainfall for Montreal. With the GATE measurements for rain 


rates associated with satellite-derived cloud areas and the 
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model output, rainfall rates were estimated by multiplying 
the two values. The most accurate estimations were for the 
cumulus clouds in the warm air Rasses OGGULE Led in gune, the 
cases the model was designed to handle. Wylie concluded 
that in order to estimate rainfall in all geographical areas 
and seasons a more sophisticated model would be needed. 
Negri and Adler (1981) did one case study of fifteen 
thunderstorms in the Oklahcma, Arkansas, and Missouri area 
omeez4 April 1975. They used radar data for ground truth and 
had special 5 minute GOES-E satellite passes over the area 
of interest. They were able to determine zhat the precipi- 
tation began falling, as indicated by radar data, for cloud 
top temperatures ranging from 229K Oo u2o0K) (~48eC Loe 139C) . 


The mean cloud tcp temperature value was 247K (-26°9C). 
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A. INTRODUCTICN 

The data set assembled for this study consists of collo- 
cated GOES-E satellite data and Service-A hourly surface 
observations for the southeastern United States during 
August 1979. The GOES-E data consists of 10x 10 pixel 
matrices of visual and infrared satellite data centered over 
each of 137 surface staticns (Fig. 8) all south of YOON, 
The satellite data are measured with the Visual Infrared 
Spin Scanned Radiometer (VISSR) which have subsatellite 
point spatial resolutions of 1 and 7 km for the visual and 
infrared channels, respectively. The GOES-E navigation was 
completed by Man-computer Interactive Data Access System 
(McIDAS) at the University of Wisconsin using the full reso- 
futicn visual data, with an accuracy of 1-2 pixels (1-2 km). 
The full resolution visual data were averaged to a 7 km res- 
olution, to equal the infrared data resolution. The visual 
and infrared digital counts range from values of 0-255. The 
10 x 10 pixel GOES-E visual and infrared satellit=s data each 


cover an area 45S nmi x 45 nmi at 30°N (60 nmi x 60 nmi at 
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42°ON). The Seéervice-A hourly Beports  cetal 70,623 
observations. No Service-A specials or record-sp¢cials are 


included. 


sgusbdury 





Figure 8. G2ographical Locations of Service-A Station 
Report Data 


The data are divided into two no-precipitation catego- 


ries (Table VIII) and seven precipitation categories (Table 
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IX) to investigate precipitation specification, convective 
versus continuous precipitation specification, and qualita- 
tive specificaticn of light versus moderate/heavy precipita- 
ON . The pixel array size is also varied from the 10 x 10 
array size to an 8 x 8, a6x 6, anda 4 x 4 array size to 
investigate the differences in the data resulting from vari- 
ous resolution sizes within a particular weather condition 


@imass: fication. 


B. CATA SORT 

For the combined visual and infrared threshold specifi- 
@aeion Of precipitaticn, satellite data for 1200-2000 GMT, 
Seeresponding to 0800-1600 EDT, were sorted into precipita- 
fmen and no-precipitation groups (Fig. 9). The 0800-1600 
EDT interval was chosen te avoid distortion of the visual 
satellite data due to a low solar elevation angle. The vis- 
ual data were normalized and converted to albedos based on 
the work done by Muench and Keegan (1979). This scheme cor- 
rects for the varying zenith angle as well as adjusting the 
visual satellite data for anisotropic scattering as related 
to the zenith angle. {See Appendix A for further specific 
information ccncermning the Muench and Keegan normalization 


scheme.) 
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States, Aucust 1979, 


HOURLY OBSERVATIONS (70, 623) 


SURFACE DATA 
Fast CoAsST AND GULF CoaASsT UNITED 






SERVICE - A 









DAYLIGHT TIME CHECK 
(AVOID LOW SOLAR ANGLE FOR 
VISUAL SATELLITE DATA) 


IncLupe 1200-2000 GMT 
(0800-1600 EDT) reports 














ISOLATE APPLICABLE WEATHER 
AND CLoUpD CASES FOR STUDY 
(see TaBces VIII anp IX) 









DELETE REPORTS WITH ANY 
ZERO VALUES IN SATELLITE 
VISUAL AND INFRARED DATA 










NORMALIZE VISUAL DATA 
CONVERT VISUAL DIGITAL 
COUNTS TO ALBEDOS 
(MUENCH AND KEEGAN, 1979) 










COMPUTE MEANS AND STANDARD 
DEVIATIONS OF THE ALBEDOS 
AND CLOUD TOP TEMPERATURES 

FOR EACH 10 x 10, 8 x 8, 6 x 6, 

AND 4 x 4 PIXEL ARRAY SIZE 











CALCULATE MEAN AND STANDARD 
DEVIATION OF THE MEANS AND 
STANDARD DEVIATIONS FOR EACH 
WEATHER/CLOUD CLASSIFICATION 
AND PIXEL ARRAY SIZE 






Figure 9. Flew Chart of Data Processing 
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COLLOCATED SATELLITE DATA 


VrisuaL 10 x 10 INFRARED 10 x 10 
PIXEL SIZE ‘PIXEL SIZE 





While alkedo values cannot exceed 1.00, the Muénch and 
Keegan (1979) scheme allows the values to overshoot 1.00, up 
+9 a value of 1.20. Therefore, the visual satellite values 
are not true albedos, but estimated albéedos. The extended 
visual normalized data scale was used to facilitate compari- 
son of the results in this effort to the most extensive bi- 
spectral threshold precipitation specification of Muench and 
Keegan (1979). The Muench and Keegan (1979) normalization 
scheme specifies that any computed albedo greater than 1.20 
be set equal to 1.20 to limit the unreasonably large values. 
Similiarly, the scheme specifies computed albedos less than 
0.15 be interpreted as the ground or water surface reflec- 
tance and the value 0.00 be assigned. The infrared data 
were processed in digital counts and converted to cloud top 
Semperatures pricr to statistical computations and graphical 
displays. 

The no-precipitation data (Table VIII) are comprised of 
ome digital visual and infrared 10 x 10 pixel arrays of 
those Service-A staticn reports not showing any "R" in the 
current weather group. DiS cle NO-=DEScipltatszon group 
includes stations reporting drizzle (weather codes L-, lL, 


and Lt). 
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TABLE VIit 


Classification of Ne-Precipitation Data Groups 


Service-A 
Current Weather/ 


Category Name Cloud Group Reports 
no eR 
2A No-Precipitation, Cloud Beas 1976 
Overcast Ceiling 300, 030, 003, 
S076 230, 013, 
ws, 103s, 203, 
Wis, 25, 213, 
223 
no wR tT 
2B No-Precipitation, Cloud Groups 7350 
Overcast and Broken (above eae) 
and 200, 02 
Wipe 20, 002, 
OZ e202, V2, 
W222 go 2 22 


The no-precipitation data are divided into two catego- 
ries, overcast ceiling (category 2A) and overcast and broken 
ceiling (categcry 2B). Cloud cover is based on the three 
feg2t cloud group in the Service-A surface observation. The 
first digit indicates the amount of low clouds, where 0 is 
defined as clear, 1 is scattered (one-eighth to four- 
@ighths cloud cover), 2 is broken (five-eighths to seven- 
eighths cloud cover), and 3 is overcast (¢ight-eighths cloud 
cover). The SSeend andeenird “drig@treandncacts the amount of 
middle and high clouds, respectively. The same 0-3 values 
defined fer the low clouds are used for middle and high 


cloud amount. 
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The precipitation data (Table IX) are comprised of the 
visual and infrared 10 x 10 pixel arrays of those Service-A 
Seaticn reports showing any “R" ain the current weather 
group. Two precipitation observations were excluded fron 


+he data set because each report also indicated clear skies. 


TABLE IX 
Classification of Precipitation Data Groups 


Siervise—-A 
Category Nane Current Weather Reports 
1 Precipitation any, A“ 5 38 
1A Precipitation RY, She and -Overncast 329 
Overcast Ceiling oa eae Kas defined in 
Table VIII category 2A) 
1B Precipitation anvouk! cand Overcast 534 
Overcast and and broken ee Lage as 
Breken Ceiling defined in Table VI 


category 2B) 


aC Continuous hoe hg Re WZ 
Precipitation 

1D Convective RW-, RW, RW+,TRW-, 426 
EReclpae dic lon ine Wee ee LX Rook 

1E Light — R=, ee. CRW, TR= 464 
Precipitation 

1F Moderate/H2avy R, Rt, RW, RW+, TRW, 74 
Precipitation PRN, TR, TR 


The general precipitation data (category 1) are divided 


mmco SiX groups: precipitation overcast ceiling (category 


48 





1A), precipitation overcast and broken ceiling (category 
ie, continuous (category 1C), convective (category 1D), 


light (category 1E), and moderate/heavy (category 1F) pre- 


eepitation. These seven precipitation groups ar= used to 
investigate precipitation specification, convective versus 
continuous precipitation specification, and qualitative 


specification of light versus moderate/heavy precipitation. 


Ce SrALLSTICAL TREATMENT 

The means and standard deviations of albedos and cloud 
top temperatures of each 10 x 10 pixel array for the weather 
types listed in Tables VIII and IX were calculated. Means 
and standard deviations of albedos and cloud top tempera- 
tures were also calculated for the 8 x 8, 6x 6, and 4 x 4 
pixel arrays centered over the surface station. Theme, x0, 
6 x 6, and 4 x 4 pixel arrays are equal to 36 Rmi x 36 nmi, 
27 nmi x 27 nmi, and 22 nmi x 22 nmi at 309°N respectively. 
Variation of the digital satellite areal coverage is used to 
investigate the differences in the statistics due +o the 
chosen resolution size. 

The data sets in Tables VIII and IX are represented, 
first, by the mean and standard deviation of the resolution 


cell means and standard deviations. Second, these data sets 
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are represented by the distributions of the mean cloud top 
temperatures and albedos where the mean cloud top téempera- 
tures are sorted into ten Kelvin (K) intervals and the mean 
albedos are sorted into 0.10 intervals. These representa- 
tive statistics and distributions are calculated for the 
four pixel array sizes. 

Micmsedatistical “and distribution tesults for differing 
resolution sizes, bi-spectral threshold specification of 
precipitation, and separation of light from noderats/heavy 


precipitation are discussed in Chapter IV. 
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A. INTRODUCTION 

The figures presented in this chapter display the dis- 
tributions of the grand means of the resolution cell means 
of albedos and cloud top temperatures for the data sets 
listed in Tables VIII and IX for the four array sizes. The 
mean values are sorted into ten Kelvin intervals and 0.10 


estimated albedo intervals. 


B. RESOLUTICN 

The effect of satellite resolution in répresenting gen- 
eral precipitation (category 1) and no-precipitation cver- 
cast (category 2A) data are explored for four resolution 
sizes. ThewECUL SilZes are 10% 10, 8 x 8, 6 x 6, and 4 x4 
and ar¢ approximately equal to areas of 2025 nmi?, 1296 
nmi2, 729 nmi@, and 484 ni? at 309N respectively. 

The general precipitation (category 1) and no-precipita- 
tion overcast (category 2A) data were chosen for study 
because, while they represent two different weather condi- 
tions, their albedo and cloud top temp2rature distributions 


have the largest amcunt cf overlap when compared to any 
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Seeer Palit Of precipitation versus no-precipitation data 
sets. The possibility arises that statistical differences 
in the four resolution sizes might be sufficient or comple- 
ment other informtion in delineating thes? two weather 
@enditions. 
1. Precipitation Data 

The general precipitation (category 1) data are con- 
prised of 329 overcast ceiling reports (61%), 205 broken 
ceiling reports (38%), and & scattered ceiling reports (1%). 

a. Mean Statistics 

The precipitation data (category 1) differences 

between the means of the cell means visual and infrared 10 x 
10 and & x 4 array sizes are 0.035 and 2.2K, respectively 
(Table X). The trend of the mean of the means is toward 
higher alkedo values and cclder cloud top temperatures with 
the decreasing area or array size. The standard deviations 
of the means simliarly show an increase in the albedo, 
mero, and cloud top temperature, 1.0K, - from the 10 x 10 
array size to the 4 x 4 array size. 

b. Standard Deviation Statistics 

The standard deviation statistics display the 


Opposite trend with decreasing area as the mean statistics. 


Se 





TABLE X 


Precipitation Data Statistics for Four Array Sizes 


Overcast, Broken, and Scattered Ceilings 


10 x 10 8 x 8 6 x 6 4x4 
Mean of (VIS) 579 oo | 2603 ~614 
Means 

(IR) 253.0K Dee 25 We tk 20h o Kk 

(~ 20°C) (-2 19°C) (-22°C) (~22°9C) 
Standard (VIS) o 211 2 214 0219 weer 
Deviations 
of Means (IR) 20.7K Ze OK lee SK oN. TK 
Mean of (VIS) eS - 161 2144 ~124 
Standard 
Deviations (IR) 10.4K 9.1K TOK DOK 
Standard 
Deviation (VIS) - 078 078 5075 O72 
of the 
Standard 
Deviations (IR) ero 6.7K 6.0K 4.9K 


The means and standard deviations of th2 standard deviations 
decrease in the visual and infrared values with decreasing 
area (Table X). The differences between th2 10 x 10 and 4 x 
4 array sizes visual and infrared means of the standard 
deviations are 0.049 and 4.8K, respectively, and the stan- 
dard deviaticns of the standard deviations are 0.006 and 


weik, cespectively. 
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c. Distribution Discussion 

The distributicns of the precipitation data are 
Shown in Figs. 10, 11, 12, and 13. There is a discernible 
shift toward higher albedcs and colder cloud top t¢empera- 
tures of the 2% and 3% frequency isopleth with decreasing 
acray size. This upward shift is also reflected in the mean 
of the means (Table X). The appearance of the 5% frequency 
tsopleth in the 6x 6 and 4 x4 array sizes at high albedos 


and cold cloud tcp temperatures highlights the shift. 
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310 300 290 280 270 260 250 240 230 220 210 # 200 


ee al CLOUD TOP TEMPERATURE (K) 


Figure 10. Precipitation Data for 10 x 10 Array Size (The 
méan, .597/9 and 253.0K, interval is boxed. The 
2% and 3% frequencies are for 11 and 
occurrences, respectively.) 
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ALBEDO 





310 300 290 280 270 260 250 240 230 220 6 ©210 8 200 


CLOUD TOP TEMPERATURE (K) 


meagre 11. Precipitation Data for 8 x 8 Array Size (The 
mean, .991 and 252.2K, interval is boxed. The 
2% and 3% frequencies are for 11 and 16 
occurrences, respectively.) 
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310 300 290 280 8 270 260 250 240 230 220 210 #8200 
CLOUD TOP TEMPERATURE (K) 
Figure 12. Precipitation Data for 6 x § Array Size (The 
Bor ee and 251.4K, interval 1S boxed. The 
¢ 


Lk, and 5% frequencies are for 11, 16, and 
27 occurrences, respectively.) 
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ALBEDO 





310 300 290 280 270 260 250 240 230 220 8 210 8 200 


CLOUD TOP TEMPERATURE (K) 


Figure 13. Precipitation Data for 4 x 4 Array Size (The 
mean, .614 and 250.8K, interval is boxed. The 
2m, 5%, and 5% frequencies are for 11, 16, and 
2) cccurrences, respectively.) 


The distributions of the 10 x 10 and 8 x 8 array 
sizes (Figs. 10 and 11) are unimodal while the 6 x 6 and & x 
4 array sizes (Figs. 12 and 13) appear to b2 more bimodal. 
The four array sizes were tested for a Gaussian distribution 
With the Chi-square test and all failed az any confidence 
level.: Therefore, differences in the four resolutions can- 
not be adequately tested by well defined statistical methods 
based on an assumed normal distribution. 

The similiarities between tha 10 x 10 and 8 x 8 


array sizes (Figs. 10and 11) and the 6 x 6 and 4 x 4& array 
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sizes (Figs. 12 and 13) are further illustrated in Fig. 14. 
A diagonal cut is plotted for each of the four array sizes 
where the lines plotted are shown as a dashed boxed area in 
Figs. 10-13. The diagonal cut reveals the close agreement 
between the 6 x 6 and 4 x 4 array sizes along the line. The 
10 x 10 and 8 x 8 array size lines follow the same general 
«rend but do not cincide as closely as the 6 x 6 and 4 x 4 
array size lines. 

The chosen diagonal line results in the 8 x 8 
array size distributicn appearing more smoothed than the 10 
x 10 (Fig. 14), as there is no relative minima at the 4.0 
interval for the 8 x 8 array Size. The 10 x 10 array size, 
with the greater areal extent and therefor= more averaging 
of differing clouds and clear areas, is expected to possess 
the "smoothest" appearance, the lowest number of relative 
maxima and minima of the four array siz2s. However, the 8 x 
8 array Size actually displays the fewest relative maxima 
and minima along the chosen diagonal line. The smoother 8 x 
8 array size cannot be explained in terms of significant 
differences in the number of cases of different ceiling 
types or different precipitation types occurring in interval 


4.0 between the four array sizes. Quite simply, the 


7. 


smoother 8 x 8 array size 


apparently results from the sort- 
ing intervals chosen for the distributions. 
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Figure 14. Precipitation Array Size Distributions Along 
Diagonal Line (The line represents the 10 x _10, 
dotted line tke 8 x 8, dashed line the 6 x 6, 
and dash dot line the 4 x 4 array size.) 
The elongated shape of all four precipitation 
distributicn array sizes (Figs. 10, 11, 12, and 13) reveal 
the variaticn in the areal amount of cloudiness and precipi- 
mares CN. 


The distributions range from high albedos and cold 
cloud top 


temperatures (indicative 


of satellite fields of 
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view filled with precipitating clouds) to low albedos and 
warm cloud top temperatures (indicative of satellite fields 
of view partially filled with precipitating clouds). The 
visual and infrared satellite data distributions in Figs. 
10, 11, 12, and 13 agree with the elongated shapes of Platt 
(1981) for his cloud classifications of cumulus, frontal, 
and jetstream cirrus clcuds and agres with Coakley and 
Bretherton (1982) for their general clouds present in a 1000 
km@ Pacific Ocean area. 
d. Summary 

A satellite field of view filled with a precipi- 
tating cloud is expected to have high2r albedo and colder 
cloud top temperature values than a partially filled field 
of view. Additionally, the filled field of view would have 
@ more uniform texture, as reflected in variance or standard 
deviation values, RinameaePacctalive filled field of view. 
The statistics discussed in this study confirm these expec- 
meerens for this data set. As the array size or fieid of 
view is decreased, the mean statistics increase while the 
Standard deviation statistics decrease (Table X). 

Fer the precipitation data (category 1), there 


are significant differences between the four resolution 
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sizes. These differences are reflected in the upward trend 

= 
in albedos and cclder cloud top temperatures of the mean and 
the standard deviation of the means with decreasing area or 
array size (Table X). The reverse trend is found in the 
mean and standard deviation of the standard deviations. 

Tke distributicns have significant differences 
also. The relatively ccarse resolution 10 x 10 and 8 x 8 
array sizes have a unimodal distribution while the rela- 
tively fine resclution 6x 6 and 4 x 4 have a bimodal 
mestribution. 

The four array sizes discussed vary in their 
statistics and distributicns in representing the precipita- 
fen (Category 1) data. The choice of satellite resolution 
for representation of the precipitation data will influence 
comparison of these data with other dat3. Therefore, for 
the remainder of this study the precipitation data for all 
Classifications will be discussed using both the 10 x 10 and 
4 x 4 array size. 

One additional topic to explore is that a number 
of cases with albedos less than 0.40 appear in the precipi- 


tation data (category 1) in all four size distributions 


(Figs. 10, 11, 12, and 13). These low albedo values suggest 
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the possibility that there might be a consistent low bias in 
the normalization scheme. Butea closer look at the 
individual reports with albedos less than 0.40 show no 
pattern involved with either the GMT hour or the longitude 
or latitude of these staticn reports. Further analysis of 
these low albedo precipitation reports are discussed in the 


light precipitation secticn (IV.E.). 


a. Mean Statistics 

For the no-preéecipitation overcast cases (cat- 
egory 2A), differences between the means of the means visual 
and infrared 10 x 10 and 4&4 x 4 array sizes are 0.011 and 0.3 
K, respectively (Table XI). The 0.3K infrared difference is 
meaan the 9.5K Roiseualenaa of the VISSR infrared sensor. 
The trend of the visual mean of the means is upward with the 
decreasing array size. The standard deviations of the means 
show an increase in the altedo, 0.014, and cloud top temper- 
ature, 0.8K, from the 10 x 10 to 4 x 4 array sizes. Once 
again, both mean statistics have an upward trend with 


decreasing array size. 
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TABLE XI 


Ne-Precipitation Data Statistics for Pour Array Sizes 


Mean of (VIS) 
Means 
(IR) 


Standard (VIS) 
Deviations 
of Means (IR) 


Mean of (VIS) 
Standard 
Deviations (IR) 


Standard 
Deviation (VIS) 
of the 
Standard 
Deviations (IR) 


b. Standard Deviation Statistics 


Overcast Ceiling 


eee 2 
50 4K 


- 056 


4.8K 


Cenversely, 


have a downward trend with decreasing area size. 


ferences between 


mmacerared 


2.5K, respectively. 
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means of the standard 


The differences in 
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Statistics 


The dif- 


Visual and 


are 0.035 and 


the standard devia- 


tions of the standard deviations are 0.007 in the visual and 


1.7K in the infrared 


values. 
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c. Distribution Discussion 

The distributions of the no-precipitation over- 
cast data are shown in Figs. 15, 16, 17, and 18. These dis- 
tributions are quite different from the precipitation 
Meeetcibutions (Figs. 10, 11, 12, and 13). As expected, the 
no-~precipitaticn overcast data are clustered at the low 
albedo and warm cloud top temperature values. The 5% fre- 
quency isopleths in Figs. 15, 16, and 17 show a grouping of 
the data at albedos ranging from 0.30 to 0.50 and cloud top 
temperatures frem 280K tc 290K. A bimodal distribution 
appears in the finer resolution 8 x 8, 6x 6, and 4 x 4 
array sizes (Figs. 16, 17, and 18). This shifting of the 
no-precipitaticn overcast data into two rslative maxima for 
the thre2 smallest array sizes is the sole significant dif- 
ference in the four distributions. 

d. Summary 

The no-precipitation overcast data display an 
upward trend in the two mean statistics with decreasing 
array size, although the 0.3K mean of the means infrared 
Maeererence is not significant. Conversely the two standard 
deviation statistics decrease with decreasing array size. 


These trends are consistent with the expected statistical 
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ALBEDO 


Figure 15. 


ALBEDO 


moegure 16. 





310 300 290 280 270 260 250 240 230 220 210 #200 
CLOUD TOP TEMPERATURE (K) 


Ne-Precipitation Overcast Data for 10 x 10 Array 
Size (The mean, .407 and 272.6K, interval is 
boxed. The 2%, 3%, 5% and 7% frequencies are 
mom WO, oo, 99, and 138 occurrences, 
respectively.) 





310 300 290 280 270 260 250 240 230 220 # 210 = #200 


CLOUD TOP TEMPERATURE (K) 


NOsmece patatdon Overcast Daca £or 39 x 8 Array 
Size (@he mean, .411_ and 272.5K, intezval is 
boxed. The 2%, ft 5% and 7% frequencies are 


for 40, 59, 9Y, and 138 occurrences, 
respectively.) 
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ALBEDO 


meaure 17. 


ALBEDO 


Figure 18. 





310 300 290 280 270 260 250 240 =230 220 210 += # 200 
CLOUD TOP TEMPERATURE (K) 


No-Precipitation Overcast Da-a for 6 x 6 Array 

Size (he mean, .415 and 272.4K, interval is 

boxed. The 2%, 3%, and 5% frequencies are for 
0, 59, and 99 occurrences, respectively.) 
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310 300 290 280 270 260 250 240 230 220 ~=© 210 200 


CLOUD TOP TEMPERATURE (K) a 


No-Precipitation Overcast Data for 4 x 4 Array 
size “ihe mean, -418 and 272.3K, interval is 
boxed. The 2%, 3%, and 5% frequencies are for 
4Q, 59, and 99 occurrences, respectively.) 
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trends discussed in the precipitation section (IV.B.1.). 
However, in these data, Bass is a greater similiarity in 
the statistics for the four sizes because the ceilings are 
all cvercast reports. 

The distributions (Figs. 15, 16, 17, and 18) 
provide visual confirmaticn of the similiarities between 
each of the four array sizes. With the exception of the 
second relative maxima at albedos of 0.30 to 0.40 and cloud 
top temperatures of 260K tc 270K appearing in the 8 x 8, 6 x 
6, and 4 x 4 array sizes (Figs. V6rmiv> and is), the four 
distributions are nearly identical. Because there are sta- 
tistical differences and distributional differences in the 
four array sizes, the 10 x 10 and 4 x 4 arrays sizes will be 
used to represent the two no-precipitation data categories. 

3. Precipitation and No~precipitation Comparison 

The possibility arises that statistical differences 
in the four resolution sizes might be sufficient, OE ae 
least complement other information, in delineating the pre- 
Cipitation frem the ne-~precipitation weather condition. The 
question then becomes, is there a statistic associated with 
Variation of the array size within Tables X and XI which 
differentiates the precipitation reports from the no-precip- 


itation overcast reports? 
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iPemeene Satellite data processor can vary the 
resolution size, as was dcne in this study by simply averag- 
ing different pixel array sizes, a trend in the visual and 
infrared data might be used to differentiate these two 
weather conditions. The most Significant trend difference 
in the precipitation data (Table X) ade 1o-pre Grpitacaon 
overcast data (Table XI) occurs in the mean of the means. 
Recall that the mean of the means precipitation visual dif- 
ference between the 10 x 10 and 4 x 4 array sizes was 0.035 
while the no-précipitation overcast visual difference was 
wero l 1. Similiarly the infrared differences were 2.2K and 
Meh £LOr the precipitation and no-precipitation overcast, 
respectively. The precipitation data show a greater upward 
trend toward higher albedeos and colder cloud top tempera- 
tures than the no-precipitation overcast with the finer sat- 
ellite resolution. 

Differentiation between these two data sets based on 
a comparisen cf che trend in the mean of the means are sug- 
gested by Tables X and XI. It must be emphasized that these 
tables are based on many reports and therefore reflect the 
most typical values. PGi VvVadvel ETeports within a given 


interval should be studied to provide conclusive evidence as 
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to whether these statistics can be used on a few reports to 
differentiate precipitaticn from no-pracipitation overcast 


reports. 


Ce PRECIPITATICN SPECIFICATION 

An essential difference between this specification study 
and most of those in the literature is that the distribu- 
tions of the precipitaticn and no-precipitation data sets 
are examined in detail to extract information about the 
probability of correct classifications. Only Lovejoy and 
Austin (1979) present their data distributions. The bi- 
spectral and life history method thresholds (Tables I and 
7 refer tc the typical or most common threshold values 
for precipitation, which is assumed to bé equivalent to «he 
mean of the méans in this study. Therefore, the mean of the 
means can be compared to the threshold values in Tables I 
and VII. Additionally, a bi-spectral threshold can be pro- 
posed based on the distributions and wit these distribu- 
tions the amount of overlap, or the percentage of correctly 
classified precipitation or no-precipitation ee can be 


calculated. 


68 





1. Overcast Ceilings 
Are the precipitation overcast (category 1A) and 
no-precipitation overcast (category 2A) data sets suffi- 
ciently separated to allow differentiation of th? two popu- 
lations? If so, how much overlap is there between the two 
data sets? 
a. Mean of the Means 

The mean of the means statistics for the precip- 
itation overcast versus no-precipitation overcast data, 
Table XII, show there is a .242 and 24.4K difference between 
the two populations for the 10 x 10 array size and a .254 
and 25.9K difference for the 4 x 4 array size. The respec- 
“ive differences are greater than one standard deviation of 
the means of either of the two populations. 

If the two precipitation overcast array sizes 
are compared to the no-precipitation overcast array sizes, 
the stronger trend in the mean of the means (Table XII) is 
seen in the precipitation cvercast data. While the precipi- 
tation visual and infrared values vary by 0.023 and 1.8K 
between the 10 x 10 and 4 x 4 array sizes, the no-~precipita- 
tion visual and infrared values vary by 0.011 and 0.3K. 


Both the mean c£ the means and their trends can be used to 
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TABLEX LE 


Precipitation Specification Overcast Ceilings 


bigcel.p.~ Nowereer~ Prect pl | NO-~Prec.— 
faelon Picate on eat lon plete t2on 
10x10 10 x 10 4x4 4x4 
Mean of (VIS) - 649 - 407 sor 2 ~418 
Means 
(IR) 248.2K 2 20K 246.4K LIES, 
(- 25°C) (- 19C) (-27°C) (-19C) 
Standard (VIS) - 186 20 22.01 2219 
Deviations 
of Means (I R) 18.9K 17. 8K 19.8K 18.6K 
Mean of (VIS) me 2 alee ~ 108 - 087 
Standard 
Deviations (IR) 8.9K 5.4K 4.4K Ze Ke 
Standard 
Deviation (VIS) - 078 « 056 069 049 
of the 
Standard 
Deviations (IR) O25 K 4.8K 4.0K Seti 


differentiate precipitation from no-précipitation for a 
large number of reports ina region similiar to this summer- 
time convective shower dominated area. 
b. Meéean of the Standard Deviations 

The visual andinfrared means of the standard 
deviations vary by 0.030 and 3.5K for the 10 x 10 array size 
ema by 90.021 and 1.5K forthe 4 x 4 array size. In the 
infrared values, the no-precipitation mean of the standard 


deviations have a magnitude 60% of the precipitation values 
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and produce relatively large differences. In the visual 
values, the no-precipitaticn mean of the standard deviations 
have a magnitude 80% of the precipitation values. The rela- 
tively large differences (cnly the mean of the means have a 
larger difference) suggest the use of this statistic to dif- 
ferentiate precipitation overcast EEOteeene- plreeraps+a tion 
overcast. 
c. Standard Deviation of the Means 

The differences in the standard deviations of 
the means (Table XII) are nearly equal when comparing the 
two 10 x 10 array sizes and the two 4 x 4 array sizes. The 
Visual differences are 0.019 for the 10 x 10 and 0.018 for 
the 4 x 4 array size. Similiarly, the infrared differences 
meen l.1K for tke 10 x 10 and 1.2K for the 4 x 4 array size. 
The differences in the 4 x 4 array size (0.018 and 1.2K) are 
comparable to the relatively significant differences in the 
means of the standard deviations (.021 and 1.5K). However, 
the differences are not comparable in the 10 x 10 array 
Size. 

d. Standard Deviation of the Standard Deviations 
The differences in the standard deviations of 


the standard deviations (Table XII) for the two array sizes 
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are approximately equal also. The visual differences are 
0.022 for the 10x 10 and 0.020 for the 4 x 4 array size. 
The infrared differences are 1.5K for the 10 x 10 and 0.9K 
for the 4 x 4 array size. Once again the differences in the 
4 x 4 array size (0.020 and 0.9K) are approximately equal to 
the differences in the means of the standard deviations 
(0.021 and 1.5K), particularly in the visual value. 
e. Distribution Discussion 

The distributions for the 10x 10 array size 
precipitation overcast (Fig. 19) and no-precipitation over- 
cast (Fig. 15) and the 4 x 4 array sizé precipitation over- 
cast (Fig. 20) and no-precipitation overcast ((Fig. 18) 
allow visual confirmation of the degree of separation of 
these two populations. These figures verify the separation 
between the the cccurrence maxima of the two populations 
while showing that there is overlap of some of the values in 
the two porfulations. 

The appearance of a bimodal distribution in the 
relatively fine resoluticn4 x 4 array size precipitation 
Overcast (category 1A) and no-precipitation overcast (cat- 
egory 2A) in Figs. 20 and 18 cannot be explained in terms of 


the precipitation categories (Table IX) or ceiling cover. 
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Bequre 19. 


ALBEDO 


meagure 20. 





310 300 290 280 270 260 250 240 230 220 210 = #8200 
CLOUD TOP TEMPERATURE (K) 


Precipitation Overcast Data for the 10 x 10 
Array Size (The mean, .649 and 248.2K, interval 
is boxed. he 2%, 34, and 5% frequencies are 


for 7, 10, and 16 occurrences, respectively.) 





310 =300 290 280 270 260 250 240 230 220 210 #200 


CLOUD TOP TEMPERATURE (K) 


Precipitation Overcast Data for 4 x 4 Array Size 
(Tke mean, «672 and 246.4K, interval is boxed. 
The 2% 3%, 5% and 7% frequenciss are for 7, 10, 
16, and 23 occurrences, respectively.) 


73 





An alternate explanation might be that these two maxima 
reflect different synoptic signatures. Four weak frontal 
systems impact this data set region during August 1979 and 
cause a change in the clceudiness and precipitation pattern 
which is normally produced by daytime heating. Recall that 
the Lovejoy and Austin (1979) data set also showed bimodal 
distributions for the cumulus no-rain reports. The bimedal 
distributions in this study may not b@? due +o a synoptic 
signature. Nonetheless, Geis  poessability should be 
investigated. 
f. Precipitation Probabilities 

PEecipleatzon EEeoObabilities (Figs. 21 and 22) 
were computed from the precipitation overcast (category 1A) 
and the no-preéecipitation cvercast (category 2A) data for the 
10 x 10 (Figs. 19 and 15) and the 4 x 4 (Pigs. 20 and 18) 
array sizes. Estimated albedes greater than 1.00 were not 
included in these prcebabilities as they accounted for only 
two and three no-precipitation overcast reports and four and 
five precipitation overcast reports in the 10 x 10 size and 
4 x 4 array size, respectively. The two probability figures 
maeacate that the 50%’ prokability cf precipitation is not a 
Simple function of mean albedo and mean ELeUud = etoD 


temperature. 
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occurs at high albedos (0.80-1. 00) aia dencord Cloud top 
temperatures (210K-23 0K) occurs at relatively low albedos 
(0.40-0.60). One exception occurs in the 0.60-0.80 albedo 
and 210K-230K cloud top temperature interval and represents 
26 precipitaticn reports cf 64 total reports. The 100% pre- 
cipitation probability at 0.00-0.20 albedo and 230K=-250K 
results from two precipitation reports. Two reports ina 
0.20 albede and 20K cloud top temperature interval aré not a 
sufficient number of reports to be a significant indication 
Sena high prokability of precipitation. The 50% precipita- 
tion probability line (Fig. 22) in the 4 x 4Y array size data 
shows the same general trend of decreasing albedo with 
decreasing cloud top temperatures. Eetias Lane soso lution 
data (Fig. 22), there is an upturning of the 50% probability 
line at the coldest cloud top temperatures, 210K-230K. 
There are 73 no-precipitaticn overcast and 80 precipitation 
overcast reports in the 4 x 4Y array size data in the 
210K-230K interval so the upturning is not the result of 
lack of data. The appearance, once again, of a greater than 
50% probability cf precipitation at low albédos, 0.00-0.20, 
between 210K-250K results from a total of four reports (one 


precipitation report of two total reports in the 230K-250K 
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interval and two precipitation reports of two reports in the 
210K-230K) and is not a significant indication of high pre- 
@mpitation prokability. 

If a straight line is drawn to represent the 50% 
probability line, the lines for the two array sizss are 
nearly coincident until the 230K cloud top temperature is 
reached. For the purposes of this study, a simple linear 
function bi-spectral precipitation threshold based on a 50% 
probability of precipitaticn can be approximately defined as 
extending from 1.00 albedo and 290K cloud top temperature to 
0.60 albede and 210K cloud top temperature. 

Ccmparison of this linear bi-spectral threshold 
With the Muench and Keegan (1979) threshold (Fig. 2), shows 
the proposed threshold has an albedo approximately 0.10 
smaller at corresponding cloud top temperatures. The 50% 
probability asymptote at the warm cloud top temperatures 
shown in the Muench and Keegan (1979) results (Fig. 2), are 
not shown in Figs. 21 and 22 due to lack cf reports in these 
values. The lower albedo values for this linear bi-spectral 
threshold may te associated with the dominance of convective 
precipitation (426 of 538 reports) in the precipitation 
data. The effect of convective precipitation reports on the 


satellite albedo values is discussed in Section IV.D. 
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g. Summary 

Phe wewoesidntitcane Statistics for precipitation 
specification of overcast ceiling reports for both array 
sizes are th2 mean of the means and the mean of the standard 
deviations (Table XII). The probabilites in Figs. 21 and 22 
make use of the mean of the means only. Inclusion of the 
mean of the standard deviations in the precipitation and 
no-precipitaticn distribution flots may more distinctly 
define the two weather data types. The question is, how to 
graphically display four variables (i.e. four dimensions) in 
one plot? 

One solution is to find a three-dimensional plot 
that involves the four variables. The most straight forward 
approach is tc define the 50% probability for the mean of 
the means (Figs. 21 and 22) in terms of a surface. A plane 
would be the simplest surface choice. The equation of a 
line perpendicular to the 50% probability plane intersecting 
its midpoint in the distribution planes shown would then be 
calculated. All of the precipitation and no-precipitation 
points would ke projected onto the line and the line would 
become the x-axis in a new plot. Thus this x-axis reflects 


*he visual and infrared mean of the means. The visual and 
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infrared mean of the standard deviations would define the 
y-axis and z-axis, respectively. 

The defined pict for the precipitation and no- 
precipitation data would show the relative dependence of the 
data on the four variables. If the plots produced distinct 
groups for the precipitation and no-precipitation data, ne¢w 
precipitation frobabilities would be calculated. The varia- 
tion with satellite resclution size of the precipitation 
probabilities would then have to be reconsidered with the 
new data. The data processing described is beyond the scope 
of this particular research effort and 1s recommmended for 
future investigation. 

PEccupreatmoneopecit 1Gation for overeast  ceii- 
ings can te delineated by the values for the meen of the 
means and the mean of the standard deviations for any array 
size discussed in this study. The mean of the mean values 
when used with Figs. 21 and 22 will indicate the probability 
of precipitation, given a similiar time of year and climato- 
logical area. A simple linear bi-spectral threshold, based 
on a 50% probability of precipitation, is defined approxi- 
mately as extending from a 1.00 albedo and a 290K cloud top 
temperature to a 0.60 albedo and a ZACK “cloud -t£op 


temperature. 
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Ceilings 
Once again the question arises, are the precipita- 
¢ion overcast and broken (category 1B) and no-precipitation 
overcast and kroken (category 2B) data sets sufficiently 
separated to allow differentiaticn of ths two populations? 
If sc, how much overlap is there between the two data sets? 
a. Mean of the Means 

The mean of the means statistics for the precip- 
itation overcast and broken (category 1B) versus no-precipi- 
tation overcast and broken (category 2B) data, Table XIII, 
show there is a .309 and 27.0K difference between the two 
populations for the 10 x 10 pixel size and a .339 and 29.1K 
difference for the 4 x 4 fixel size. The respective differ- 
ences are approximately ¢qual to one and one-half standard 
deviations of the means of either of th2 two populations. 

For these data, the differences in the trends of 
the mean of the means between the 10 x 10 and 4© x 4 array 
Sizes for the precipitaticn overcast and broken reports are 
more dramatic than for the precipitation overcast reports. 
There is a 0.036 visual and a 2.2K infrared difference in 
the precipitation overcast and broken data and only a 0.006 


Visual anda 0.1K in the no-precipitation overcast and 
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ES pets ek eek 


Precipitation Specification Overcast and Broken Ceilings 


Precipi- Nome rect Precipi= No~Preci- 
Gael oC} Died < LOnm cacu@n Dicat2on 
10 x 10 19 x 10 ek ee 
Mean of (VIS) - 580 ser 20016 ea 
Means 
(IR) 252. OK 279.8K 250.6K 279.7K 
(~20°C) (+ 79°C) (-23°9C) (+7°C) 
tandazrd (VIS) o2ui0 aeio ZL 2 On), 
Deviations 
of Means (IR) 20.5K 16.6K Zalver on 17.8K 
Mean of (VIS) 513 - 128 e124 094 
Seandard 
Deviations (IR) 10.4K ae Oreo Sy A 
Seandard 
Deviation (VIS) - 078 - 064 072 0 0 
of the 
Standard 
Deviations (IR) 720K eK 4.9K 3.7K 


broken data between the 10x 10 and 4 x 4 array sizes. 
Therefore, variation of the resolution size and the trend 
for the mean of the means can be used, as well as the mean 
of the means value itself, in delineating precipitation from 
no-precipitation for these overcast and broken ceiling 
reports. 
be. Mean of the Standard Deviations 
The values of the mean of the standard devia- 


tions are important. The 0.045 albedo and 4.5K cloud top 
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temperature difference between the two data sets for the 10 
x 10 array size are significant as well as the 0.030 and 
2.4K difference for the 4 x 4 array siZe. 
c. Standard Deviations of the Means 

The differences in the standard deviations of 
the means are nearly equal when comparing the two 10 x 10 
array sizes and the two 4 x 4 array sizes. The visual dif- 
ferences are 0.021 for the 10 x 10 and 0.018 for the 4 x 4 
array size. Similiarly, the infrared differences are 3.9K 
for the 10 x 10 and 3.7K for the 4 x 4 array size. The vis- 
ual differences are 40% less than the statistically signifi- 
cant visual mean of the standard deviation differences and 
the infrared differences are comparable to the infrared mean 
of the standard deviation differences. 

d. Standard Deviation of the Standard Deviations 

The differences in the standard deviations of 
the standard deviations are approximately equal. The visual 
differences are 0.014 for the 10 x 10 and 0.015 for the 4 x 
4 array size. The infrared differences are 1.3K for the 10 
x 10 and 1.2K for the 4 x YW array size. Both the visual and 
infrared differences are 50% less than the respective sig- 


nificant differences in the mean of the standard deviations. 
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e. Distribution Discussion 

Tees distLibucions fer the 10 x 10 array size 
precipitation cvercast and broken (Fig. 23) and no-precipi- 
taticn overcast and broken (Fig. 24) and the 4 x 4 array 
size precipitation overcast and broken (Fig. 25) and no-pre- 
cipitation cvercast and broken (Fig. 26) ailow visual con- 
firmation of the relatively larger degree of separation 
between these two populations compared to the overcast ceil- 
ing data. Tkese figures once again also show overlap 


between the two populations. 
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Figure 23. Precipitation Overcast and Broken Data for 10 x 
Atray (The mean, .580 and 252.8K, interval i 
boxed. The 2% and 3% frequencies are for 11 an 

16 occurrences, respectively.) 
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ALBEDO 


Figure 24. 


ALBEDO 


Pegure 25. 
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CLOUD TOP TEMPERATURE (K) 


No-Precipitaticn Overcast and Broken for 10 x 10 
pemee (The meen, rae and a oe 8K, interval is 


boxe The | h%, 5%, and 7% frequencies are 
for 147 ete 63. and’515 occurrences, 
poem le aly.) 





310 300 290 280 270 260 250 240 230 220 210 = 200 


CLOUD TOP TEMPERATURE (K) 


Precipitation Overcast and Broken Data for 4 x 4 
boxed, (The mean, .616 and 250.6K, interval is 
sees The 2%, 3%, and 5% frequencies are for 
6, and 27 occurrences, respectively.) 
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Figure 26. No-Precipitaticn Overcast and Broken for 4 x 4 
RELee on = ae eel rae gong 219e kee eAterval is 
boxe ; and 12% frequencies 
are for ay ; eq a Be "515, and 883 occurrences, 
respectively.) 


f. Precipitation Frobabilities 

Precipitation probabilities (Figs. 27 and 28) 
were computed from the precipitation overcast and broken 
(category 1B) and the no-precipitation overcast and broken 
(category 2B) cata for the 10 x 10 (Figs. 23 and 24) and the 
4 x 4 (Pigs. 25 and 26) array sizes. As in the overcast 
ceiling terorts, there is a greater than 50% probability of 
precipitation at the low albedo values, 0.00-0.20. Ens sh acge 
27, the 75% probability between 230K-250K interval results 


meom three precipitation reports of four total reports. The 


85 





overcast ceiling reports account £Or two oe the 
precipitation reports. Therefore, one precipitation and one 
no~precipitaticn broken report have been added to the inter- 
val. Similiarly in Fig. 28, the 67% probability in the 
210K~-230K interval results from two precipitation reports of 
three total referts. The overcast ceiling reports account 
Bot two of the precipitation reports. Thus, one no-preéecipi- 
tation report has been added to the interval. A few reports 
in these low albedo and cold cloud top temperature intervals 
are producing misleadingly high precipitation probabilities. 
The 50% probability line is nearly constant, at a 0.80 
albedo, at all cloud top temperatures for the overcast and 
broken ceiling data. 

The same treatment for displaying the four vari- 
ables in a three dimensional diagram suggested in the previ- 
ous section (IV.C.1.), is recommended for these data as 
well. The new displays then could be used to calculate pre- 
Cipitation probabilities that migkt allow a more accurate 
Be-spectral threshold specification of pracipitation. 

ge Summmary 
As for the overcast ceilings, precipitation 


specification of the overcast and broken ceilings can be 
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Figure 27. Precipitation Data Probability 10 x 10 Array 
Size (Overcast and Broken Ceilings) 
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Mmeoure 28. Precipitation Data Probability 4 x 4 Array Size 
(Cvercast and Eroken Ceilings) 


delineated by the values for the mean of the means and the 


mean of the standard deviations for any array size discussed 
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me this st udty. The mean of the means values when used with 
Figs. 27 and 28 will indicate the PrObanpr lity OL precipita=- 
tion, given a similiar time of year and climatological area. 
A very simple threshold fcr these data is actually dependent 


upon the visual data value being greater than 0.80. 


D. CONVECTIVE VERSUS CONTINUOUS PRECIPITATION 

The physical processes involved in convective precipita- 
«ion, cr shewers, are different from those processes usually 
involved in continuous precipitation. The inherent differ- 
ences in the precesses might result in a statistical separa- 
icon in the thresholds between the two types Of 
precipitation. For this data set, 426 of the 538 cases are 
classified as convective by the surface weather report. 
With the data set consisting of southeastern United States 
Stations in August, the deminance of the precipitation cases 


by convective reports 1S not surprising. 


The statistics for the convective and continuous 
precipitation (Table XIV) indicate differences in the means 
eeetne cell means of .013 and 5.7K for the 10 x 10 array 
size and differences of .015 and 6.8K for the 4 x 4 array 


size. As mentioned in the life history methods, convective 
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precipitation is 


higher albedos and colder cloud top temperatures. 


data, the 


colder cloud tcp temperatures. 


bus clouds. 
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expected to be associated 


convective precipitation 


TABLE XIV 


Continuous versus Convective Precipitation Specification 
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precipitation (Table XIV) 


values, the mean of the standard deviations is th2 best sta- 


tistic for differentiation 


the 


Convec- 
tive 


of these two data 
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suggest that for visual satellite 
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for infrared values, the mean of the means is the best 
Beaclstic. The use of two different statistical types *o 
qualitatively specify continuous precipitation from convec- 
tive precipitation is unique to this pair of convective ver- 
sus continuous precipitaticn data comparison. 

The infrared mean of the cell means show that the 
convective precipitation has colder cloud top temperatures 
due to the greater vertical development of the clouds. The 
greater vertical variation of the convective precipitation 
clouds is seen in the visual "texture" or visual mean of the 
standard deviaticn statistics. This "texture" difference 
also appears in the infrared values, but it is not as sig- 
nificant as the differences in the infrared mean of the 
méans. The lack of a large difference in the visual mean cf 
the means between the convective and continuous precipita- 
tion reflects the averaging of open areas and cumulus clouds 
associated with convective precipitation. Note that for the 
fine resolution 4 x 4 array size, the visual mean of the 
means is larger for the convective than the continuous 
precipitation. 

The visual means of the standard deviations have a 


Geer erence of 0.033 for the 190 x 10 and 0.032 for the 4&4 x 4 
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array size. Contrast these difference to the visual mean of 
the mean differences of 0.013 and 0.015 for the 10 x 10 and 
4 x Yarray size, respectively. The infrared mean of the 
means have a difference of 5.7K for the 10 x 10 and 6.8K for 
the 4 x 4 array size. These differences are larger than the 
infrared mean of the standard deviations differences of 3.5K 
for the 10 x 10 and 2.2K for the 4 x 4 array size. 
2. Standard Deviation Statistics 

The standard deviation of the means and the standard 
deviation of the standard deviations display quite similiar 
and relatively insignificant differences. The visual stan- 
dard deviation of the means vary by 0.013 in the 10 x 10 and 
0.007 in the 4&4 x 4 with the infrared values varying by 1.4K 
for the 10 x 10 and 1.8K for the 4 x 4 array sizes. The 
visual and infrared 10 x 10 and 4 x 4 array size differences 
in the standard deviations of the standard deviations are 
wees, 9.017, 1.9K, and 2.0K, respectively. 

3. Distribution Discussion 

The distributions of the convective and continuous 
precipitaticn for the 10 x 10 array size (Figs. 29 and 30) 
and for the 4 x 4 array size (Figs. 31 and 32) demonstrate 


eaew Similiarity of the distributions. es. Sint il tari ty 
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disallows quantitative separation of th? convective and con- 


tinucus precipitation. 
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meguce 29. Convective Precipitation Data for 10 x 10 Array 
Size (fhe mean, .5/76 and 251.8K, interval is 
bexed. The 2%, 3%, and 5% frequencies are for 
Oem, and 2t ©CCuUrr ences, Crespectively.) 


4, Summary 





Defining the qualitative delineation of convective 
versus continucus precipitation in terms of the shift toward 
colder cloud top temperatures for the convective precipita- 
tion cases is cneée available delimiter. The second qualita- 
tive separation is in terms of the shift toward higher 


albedo values of the mean of the standard deviations for the 
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Figure 30. 
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Figure 31. 
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Continuous Precipitation Data for 10 x 10 Array 
Size (The mean, .589 and 257.5K, interval is 
boxed. The 2%, 3%, and 5% frequencies are for 
2, 3, and 6 occurrences, respectively.) 
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Convective Precipitation Data for 4 x 4 Array 
Size (The mean, .617 and 249.4K, interval is. 
boxed. The 2%, 3%, and 5% frequencies are for 
9,°13, and 21 cccurrences, respectively.) 
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megure 32. Ccntinuous Precipitation Data for & x 4 Array 
Size (The mean, .602 and 256.2K, interval is 
boxed. The 2%, 3%, and 5% frequencies are for 
2, 3, and 6 occurrences, respectively.) 


convective precipitation cases as compared to the continous 
precipitation cases. 

The convective precipitation cloud top temperature 
ma@an of 251.8K (-219C) ccmpares quite well with chose asso- 
ciated with life history methods (Table VII). The continu- 
cus precipitation cloud top temperature of 257.5K (-16°9C) 
compares well with the -16°C threshold of Wylie (1979) (a 
life history threshold for Montreal, Canada, in summertime) 
and with the -12°C bi-spectral threshold of Muench and Kee- 
gan (1979). These research EERE es focused on mid-latitude 


precipitation associated with fronts. 
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E. NAPE Sn 

As there is little statistical difference between the 
convective and continuous precipitation data, the qualita- 

\ 

tive analysis of precipitation can be based strictly on all 
light or all moderate/heavy precipitation. The light pre- 
Clipitation cases are expected to have a lower mean albedo 
and warmer cloud top temperature when compared to the moder- 


ate/heavy precipitation cases. 


1. Statistical Discussion 





The statistics fer the two precipitation intensi- 
ties, Table XV, show 0.100 and 8.6K daffarences for the 10 x 
10 array size and .117 and 8.1K differences for the 4 x 4 
array size in the means of the means. The differences are 
appreximately equal to one-half a standard deviation, an 
insufficient separation for specification of these two pre- 
Sepitation classifications. The remaining three sets of 
Statistics (Table XV) show nearly identical values in the 
light and modérate/heavy precipitation for the two array 


sizes. 
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TABLE XV 


Light versus Moderate/Heavy Precipitation Specification 


Licht Moderate/ Gabe ne Moderate/ 
Heavy Heavy 
10 x 10 19 x 10 4x4 4x4 
Mean of (WEES) - 565 - 665 2598 <I NV5 
Means 
(IR) 25 4a2K 245.6K 251.9K 243.8K 
(- 199°C) (-2 89°C) (- 21°C) (~29°9C) 
Standard (VIS) ee 0 0 \ eZ - 204 
Deviations 
of Means (IR) 202K 20.4K 21.6K 2A ik 
Mean of (VIS) se 1eZ am Be 2124 229 
Standard 
Devi eeione (2k) 10. 2K Ae tears Bio 6.1K 
Standard 
Deviation (VIS) e076 - 081 2072 -079 
of the 
Standard 
Deviations (IR) TOK TOK 4.9K recat 


2. Distribution Discussion 
The distributions of the light and moderate/heavy 
precipitation for the 10 x 10 (Figs. 33 and 34) and «he 4 x 
4 array size (Figs. 35 and 36) allow visual confirmation of 
the cverlap kEetween the two intensities. Basically with 
only 74 cases cf moderateyheavy precipitation, the frequency 
Sentours (Figs. 34 and 36) reflect a noisy distribution 


created by a sorting of just two or three reports into a 


given interval. Additicnal intensity data, specifically 


oe 
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moderate/heavy precipitation cases, are needed before 4 
threshold delineating qualitative imcensleLes. Can be 


proposed. 
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meegute 33. Light Bae cae een Data for 10 x 10 Array Size 
ee 2 ees 5965 and 254.2K, int2rval is boxed. 
The 2%, 3%) and 5 % frequencies Sie Os 9, 4, 
and 53 Beet taices: respectively.) 


For the 10 x 10 and 4 x 4 array sizes moderate/heavy 
precipitation (Figs. 34 and 36), the area encompassed by the 
7% frequency isopleth at the higher albedos and colder cloud 
top temperatures gives an indaeation Of a possible 
Separation between the twe intensities. However, the bimo- 


@al nature of the 10 x 10 distribution (Fig. 34) and the 
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Pagure 34. 
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Beagurce 35. 
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Moderate/Heavy Precipitation Data for 1 x 10 


AGEay Sa2ze At é€ m2an, .665 and 245.6K, interval 
1s boxed. he 3%, 5%, and 7% frequencies are 
for 2, 4, and 5 occurrences, respectively.) 


20 
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Page Prectpitations Data for 4 m.4 Array Size 
The mean, .598 and 251.9K, interval is boxed. 
he 2% and 3% frequencies are for 9 and 14 

occurrences, respectively.) 
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Figure 36. Moderate/Heavy Precipitation Data for 4 x 4 
Array Size (The mean, .715 and 243.8K, interval 
is boxed. he 3%, 5%, and 7% frequencies are 
BObec see, ald 5 OGCENEEeEnces, Trespeczively.) 


three relative maxima of the 4 x 4 distribution (Fig. 36) 
Suggest that there will ke a significant amount of overlap 
between these two qualitative intensities. Because the sta- 
tistics fer the two rainfall rates are nearly identical 
(Table XV), with the exception of the mean of the means, the 
method discussed in the precipitation specification section 
(I¥.C.) which involved the mean of the standard deviations 
wemeid in further differentiation of two data types cannot 
be applied to these two intensities. The delineation will 
have to be derived from the mean of the mean statistics and 


their resultant distributions. 
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3. Summary 

There are slight manifestations of the expected 
downward shift inthe albedcs and the warmer cloud top temp- 
eratures for the light precipitation relative to the noder- 
ate/fheavy precipitation. However, additional data are 
required before athreshcld can be proposed +o delineate 
fight versus moderate/heavy pieset pl tation, The 
investigation of the additional data should begin with 
consideraticn of the méan Of “the means and their 
distributicns. 

During the cell resolution discussion cof Section 
IV.B.1., a Significant number of low albedo values (less 
than 0.40 albedc) were nected. h Gemsarason of the Laght 
precipitation data (Fig. 33) and précipitation data (Fig. 
10) for the 10 x 10 array size reveals that 96% (105 of 109 
Beports) are classified as light precipitation. Furthsr 
analysis shows that 79% (86 of 109 cases) are classified @s 
meant convective precipitation. (The decomposition of con- 
vective and continuous precipitation data into light and 
modéerate/heavy are not shown.) These low albedo values are 
not unreasonartle for light convective precipitation which 


can be produced by small, isolated cumulus congestus clouds. 
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V. SUMMARY AND CONCLUSIONS 


A. DATA PROCESSING SUMMARY 

The data were first filtered for the correct time, 
1200-2000 GMT, to ensure reliable visual satellites values. 
The visual counts were converted to albedos according to the 
Muench and Keegan (1979) normalization scheme (ses Appendix 
A). This scheme corrects for Lambertian scattering as well 
as anisotropic cloud radiation. SOLtIng the GaramaccOrding 
to the Service-A reports (current weather and cloud group) 
provided seven precipitaticn groups: convective, continuous, 
light, moderate/heavy, general, overcast ceiling, and over- 
cast and t-roken ceiling precipitation and two no-precipita- 
“ion groups: no-precipitation overcast and no-precipitation 
overcast and broken. In order to investigate the impact of 
Satellite resolution on precipitation specification, the 
stated groups were subdivided into four sizes: 10 x 10 (45 x 
mnt), & & 8 (36 x 36 nmi), © x 6 (27 x 27 nmi), and 4 x & 
(22 x 22 nmi) pixel array sizes. The means of the albedos 
and cloud top temperatures for each array size for e¢éach 


weather or cloud group were calculated. The mean values 
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were sorted into 10K cloud top temperature intervais and 
0.10 estimated albedo intervals to provide a visual d¢scrip- 
Peon Of the distributions. The standard deviation of the 
means and the mean and standard deviation of the standard 
deviations for toth albede and cloud top temperature were 


also computed for each distribution. 


B. Sar ioLLCs SUMMARY 

The computed mean and standard deviation of the means 
and the mean and standard deviaticn of the standard devia- 
tions for roth albedo and cloud top temperature were evalu- 
ated for specification of precipitation, Spec2 fica. Lon. of 
continuous versus convective precipitation, and specifica- 
tion of light versus moderate/heavy precipitation. 

Pear preclpitcation Specification (Tables XII and XIII), 
the visual and infrared mean of the means produced the larg- 
est statistical differences between the precipitation and 
no-precipition data sets fer both the overcast and the over- 
cast and broken ceiling reports. For the overcast ceiling 
reports (Table XII), the separation in the visual and infra- 
red means of the means were one standard deviation. For the 
overcast and broken ceiling reports (Table XIII), the sepa- 


ration in the visual and infrared means of the means were 
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one and one-half standard deviations. The mean of the stan- 
dard deviaticns and the trend of the mean of the means also 
yielded a large difference between the precipitation and 
no-precipitation data sets. 

Fer specification of continuous versus convective pre- 
Cipitation (Takkle XIV), the infrared mean of the means and 
the visual mean of the standard deviations yielded the larg- 
est statistical differences. Tremeouvective precipl=ation 
infrared mean of the means was 6K-7K colder than the contin- 
uous precipitation value and the convective visual mean of 
*he standard deviations was 20% greater than the continuous 
precipitation value. These two statistics, however, were 
only indications of possible separation as the two data dis- 
tributions disrlayed significant overlap. 

For light versus moderate/heavy precipitation specifica- 
tion (Table XV), the mean of the means yielded the largest 
statistical differences. However, *he relatively small 
quantitative differences (one-half standard déviation) tT 
the visual and infrared values only provided an indication 
of the upward shift in albedo andthe downward shift in 
cloud top temperature for the moderate/h2avy precipitation 


relative to the light precipitation. 


10s 





CG. PRECIPITATION PROBABILITY SUMMARY 

For the overcast ceiling reports, th2 10 x 10 array sizes 
(Fig. 21) 50% probability line decreased linearly to lower 
albedos with colder cloud top temperatures, while the 4 x 4 
array size (Fig. 22) 50% probability line decreased linearly 
*o approximately 240K and then increased linearly at colder 
cloud top temperatures. For the overcast and broken ceiling 
Leports, the 10 x 10 array size (Fig. 27) 50% probability 
line was constant at 0.80 albedo while the 4 x 4 array size 
(Fig. 28) 50% probability line gradually decreased linearly 
£0 approximately 240K and then gradually increased linearly 


at cclder cloud top temperatures. 


D. DATA DISTRIBUTION SUMMARY 

The distributions of all seven precipitation classes and 
both no-precipitation classes were non-Gaussian. This 
result disagrees with the assertion of Lovejoy and Austin 
(1979) that their cumulus rain, non-cumulus rain, and non- 
Cumulus no-rain data distributions were Gaussian. , At the 
relatively fine resolution 4 x 4 array size, all seven pre- 
Clpitation classes and the no-precipitation overcast data 
displayed bimodal distributions while at the relatively 
coarse 10 x 10 array size, all displayed an unimodal 


feesceribution. 
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te CONCLUSIONS 
The conclusicns from this study of precipitation speci- 
fication are: 


e Varying the satellite data resolution from 484 nmie2 (4 
x 4 array size) Pcomcuro sesame “(10 xX 10 array size) 
results in a statistically significant difference in 
the representation of paeee Dees HOas Or NO=proci p= tation 
data. Variation in the distribution functions and the 
characteristic. means and Standard deviations with 
increasing cell size demonstrates systematic trends 
which may with further study, provide improved basis 
for rain detection and/or quantification. 


e For overcast ceiling reports, a simple jiinear bi-spec- 
tral threshold based on a 50% probability of Bea pees 
tion is defined as extending from an albedo of 1.00 and 
a cloud top temperature of 290K t9 an albedo of 0.60 
and a cloud top temperature of 210K for overcast ceil- 
ang me ports. For overcast and broken ceiling reports 
a threshold based on the albedo being greater than 0.80 
eee a 50% probability of precipitation. However, 
the preciritation probabilities varied with the satel- 
lite resclution size. 


®e The precipitation and no-precipitation data sets can be 
differentiated by the mean of the means, the mean of 
Ee standard deviations, and the trend of the mean of 

e means. 


Seo cr rerentiation of convective precipitation from con- 
e_nuOus precipitation shows pegu =: Tieolgh tne da2,fter= 
ences in the infrared mean or the means and through the 
visual mean of the standard deviations. 


e Qualitative specification of light versus moderate/ 
heavy Becca aco 2 shows some romise through a 
threshold based on the relative shift toward colder 
cloud top temperatures and higher albedos of the moder- 
Decne y Peecrpitatron Gelative *9 ‘the light precipi- 

ation. 


ie SUGGESTED FURTHER STUDY 
The recommendations for further study are: 


e The Becpes ed data processing that would display the 
visual and infrared mean and standard deviation of the 
means (IV.C.1. should be investigated for both the 
Overcast and the overcast and broken data to determine 
Teo nNemOreeciDa tation and NO-precipitation data can be 
more distinctly separated. 
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The ee eceence of bimodal distributions at the fine 
resolutions (6 x 6 and 4 x 4 array sizes) in the seven 
precipitation categories and in the no-precipitation 
overcast category should be investigated as a possible 
Synoptic Signature. 


Formal discriminant analysis should be applied to the 
relevant statistics to yield confidence levels of indi- 
vidual visual and, infrared satellite data pairs for 
PeecCiplcation specification. 


The data should be further investigatad at the 2 x 2 
eee y Size to analyze this resolution in representing 
e data. 
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AEFPENDIX A 


The Muench and Keegan (1979) normalization relates the nor- 
malized refl¢ctivity, on , to the varying solar angle and 
maximum digital counts through the reflectance term, ©, and 
ere anisotropie scattering through the Y term. Table xvI 
defines the symbols, Table xXVII lists th2 geometric identity 
equations and Table XVIII list the normalization equations. 


Fig. 37 gives an example of the normalization applied +o the 


Stated location. 
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Symbol 


ao 


rer 6A 


m4 


TABLE X Va 
List of Symbols 


Description 
GOES video count number (0-63) 


GOES video count number for perfect 
diffuse reflector and overhead sun 


Greenwich meridian time 


Distance of earth to sun 

Mean distance of earth to sun 
Julian date 

Hour angle 

Cloud reflectivity 


Arc-length observer to subsatellite 
point 


Declination of the sun 

Zenith angle of the sun 

Longitude 

Longitude of subsatellite point 
Anisotropic scattering coefficient 
Latitude 

Azimuth of the sun 


Azimuth of the satellite 
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Units 
dimensionless 
dimensionless 
hours-minutes- 
seconds 
km 
km 
dimensionless 
radians 
dimensionless 


radians 


radians 
radians 
radians 
radians 
dimensionless 
radians 
radians 


radians 





TABLE XVII 
Basic Geometric Satellite-Earth Relationships * 


Declination: 

6 = 0.408 sin [(d-81) * 27/365] 
Solar distance ratio: 

R/R, = 1 - 0.167 cos [(d-14) * 27/365] 
Hour angle: 

h = A +m - G(hours) x 1/12 
Arc-length: 

cosy = cos(A. - A) cos? 
Satellite azimuth: 

Sin ($5 - 7) = sin CA. - A)/siny 
Solar azimuth angle: 

cost = sing? siné + cos? cosé cosh 
Solar azimuth: 


Sing, = cosé sinh/sing 


= 
Angles in radians 


eA wos OX Vee 


Muench and Keegan Normalization Equations * 


i = (=) * secZ 
Cy 
46 = Ibm - 941 
= 2 
Cc, = .cos; €(G = 50)) * 1.8) 


Cy a O.7 cos ((f = 22.5) * 4) * (1 = cost) 


C, = cos” ((4$ - 70) * 1.3) 
eee. O 480, Oot (1 + cos(2*r)) +°0.20 * (C, + C,) * Cy 
T. = [1.09 - 2*(1.09 - # * x *(R/R)*)/(2 cos” = in} 


* 
Angles in degrees 


rs 





COUNTS 
20.0 Dom 30.0 35.0 40.0 45.0 50.0 55.0 


See 


10.0 


Siskin ee 
ent eee Oe. 
coe ee ee sere 


~~ 
_— 
“= es ene 


atin Oeelo. Omega) 17-0 18.0 18.0 20.0 21.0 22.0 


pescute 37. 


eee he 
Normalized Cloud Reflectivit 2H as a Function of 
Video Count and Time of Day For April 16 at 4 30N 
and 74°W, with Calibration Joouat (ce equal to 
coe (example taken from Musnch and°Keegan, 
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